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The self-affine analysis and erraticity analysis of pseudorapidity gaps are performed 
for the data of 400GeV/c pp collisions. The self-affine analysis has been shown 
to exhibit a better scaling behavior. The self-affine multifractal dimensions and 
multifractal spectrum have been obtained. The simulated results using FRITIOF 
program can not reproduce the scaling behavior. The analysis of event-to-event 
fluctuations has been performed. The increase of event-space moments C p , q (M) 
with decreasing phase-space scale is dominated by the statistical fluctuations. The 
erraticity analysis based on measuring the pseudorapidity gaps is also performed. 
The entropy-like quantities S q and E q deviate from 1 significantly, implying that 
both of them are useful to serve as effective measures of erraticity in multiparticle 
production. The In S q versus q has a quite linear behavior, but the InEg versus q 
has only an approximate linear behavior. The FRITIOF simulated results follow 
the same scaling behavior, but the deviations from the experimental data are rather 
large. 


1 The data 

In the present investigation, the angular distribution of charged particles pro¬ 
duced in pp collisions at 400GeV/c was measured by using the LEBC films of¬ 
fered by the CERN NA27 collaboration. A total of 3950 non-single-diffractive 
events (N > 4) were measured. The accuracy in pseudorapidity in the region 
of interest (-2 < 77 < 2) is of the order of 0.1 pseudorapidity units. In order 
to compare with the experimental data, we used a Monte-Carlo (MC) gen¬ 
erator FRITIOF version 7.02 and JETSET 7.3 to simulate the multiparticle 
production in 400GeV/c pp collisions. A total of 4500 non-single-diffractive 
events (N > 4) have been created. 

2 Self-affine analysis 

Since Bialas and Peschanski proposed to study nonstatistical fluctuation in 
multiparticle production by the method of factorial moments F q (S) [1], a large 
variety of experiments were performed to search for the anomalous scaling 
behavior F q (S) oc (S)~^ q (5 —► 0). The results have shown that the power- 
law behavior does not hold exactly for high energy hadron-hadron collisions 
[2]. The reason for this is that the usual procedure for calculating higher¬ 
dimensional factorial moments is to divide phase space into bins with the 
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same M in each direction. This is called self-similar analysis. However, phase 
space in high energy multiparticle production is anisotropic. Wu and Liu 
have proposed a new method to calculate the higher-dimensional factorial 
moments which is called self-affine analysis [3]. In self-affine analysis, the two- 
dimensional phase space region AijAip is divided by X V X V . The shrinking ratios 
X v and X v are characterized by a parameter H = lnA^/lnA^, called Hurst 
exponent, which can be deduced from the data by fitting two corresponding 
one-dimensional second-order factorial-moment saturation curves 

if } (Mi) = Ai - BiM~ Ci (i = ri,<p), H v „ = (1 + C v )l( 1 + C v ) (1) 

In order to obtain the Hurst exponent from NA27 data, the second order 
factorial moment for one-dimensional phase space were calculated and are 
fitted to Eq.(3), cf. Table 1. 


Table 1 The fitting parameters obtained according to (3). 


Variables 

A 

B 

C 

x7ndf 

V 

1.371 ±0.025 

0.222 ±0.018 

0.425 ±0.109 

8.449/36 

F 

1.509 ±0.018 

0.420 ±0.039 

0.057 ±0.021 

11.16/34 


From these parameter values, we obtain the Hurst exponent H vv — 0.74± 
0.07. The second order two-dimensional factorial moments are then calculated 
with the method of continuously varying scale [4] for H = 0.5, 0.74, M v = M^ 

; H = 1.0,2.0, M v = M.y H . The results are shown in Fig.l and are fitted 
both linearly (dashed lines) and quadratically (full lines). The parameters of 
quadratic fit (y(x) = a + bx + ex 2 ) are shown in Table.2. We can see that the 
coefficient b which characterizes the strength of anomalous scaling is positive 
only for H = 0.74. The coefficient c of the quadratic term, characterizing the 
degree of upward-bending is the smallest for H = 0.74. These results mean 
that the self-affine analysis exhibits a better scaling behavior. 


Table 2 The fitting parameters obtained according to (4). 


H 

a 

b 

C 

0.50 

0.181 ±0.007 

-0.0213 ± 0.0051 

0.0131 ±0.0015 

0.74 

0.144 ±0.014 

0.0243 ± 0.0075 

0.0049 ± 0.0014 

1.00 

0.194 ±0.006 

-0.0197 ±0.0031 

0.0130 ±0.0011 

2.00 

0.192 ±0.015 

-0.0062 ± 0.0020 

0.0100 ±0.0030 


In order to obtain the self-affine multifractal spectrum, the two- 
dimensional factorial moments of continuous order [5] (q from -1 to 4 with 
step 0.2) are calculated and shown in Fig.2. The dotted lines in the figure are 
from FRITIOF Monte-Carlo, the dashed lines are from another MC having the 
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same multiplicity distribution and same number of events as the experimental 
data, but no correlations. It can be seen that no intermittency behavior can 
be observed in both MC’s. 

Using the intermittency exponents <f>(q) for continuous order q obtained 
by fitting In F(q,M) versus lnM, the multifractal dimension D(q) and multi¬ 
fractal spectrum f(a) of the self-affine fractal can be calculated through the 
following relations 

t(q) = q - 1 - 0(g), D (l) = r (g)/(<? - !), a = dr(q)/dq, f(a) = qa - r(q). 

D(q) versus q is shown in Fig.4. D(q) decreases with increasing q. The self- 
affine multifractal spectrum /(a) is shown in Fig.5. It is concave downward 
with a maximum at q = 0, /(a(0)) = £>(0) = 1. These mean that multiparticle 
production at 400GeV/c pp collisions is a self-affine multifractal process. 

The index /x for Levy stable law, defined by the equation: <j>(q) / <j>{2) = 
(q 11 — q)/ (2^ — 2) can be obtained using a method proposed by Hu Yuan et 
al. [6]. From Eqs.(5) and (6) we can get the following relationship [6] 

1 — f(a) oc (B — a)^/^ -1 ), for a<B. (2) 

where B is the value of a when f(a) = 1. From Fig.4, the value of a for 
/(a) = 1 is found to be B = 1.0260. The ln(l-/(a)) versus ln(H — a) are 
shown in Fig.5. The slope C = 2.0997 ± 0.0098 and the Levy index /i = 
C/{C — 1) = 1.91 ± 0.01 are obtained through a linear fit. 

3. The analysis of event-to-event fluctuations 

The investigation of non-linear phenomena in high energy collisions has 
lasted a long time. Cao and Hwa [7]. proposed to characterize the spacial 
pattern of an event by using the horizontally normalized factorial moments. 
However, their method is meaningful only for high multiplicity events [8]. 

When the event multiplicity N is low and the number of bins is high, 
only a few events have n m > q , so the statistical fluctuation may be very 
large and very little information can be obtained. In Fig.6 the event space 
moments In C p>q (M) versus In M from NA27 data are plotted on the left side 
and In C p>q (M) versus lnC^^Af) on the right. In order to see whether the 
variation of the In C p ^ q (M) versus InM is caused by the statistical fluctua¬ 
tions, we created a MC event sample which has same multiplicity distributions 
within A?/ as experimental data, but the particles are randomly distributed 
with equal probability. The results are shown in Fig.6 as dashed lines. We 
can see that the increase of the moments C p>q (M ) with decreasing phase-space 
scale is indeed dominated by the statistical fluctuations. 
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In order to get the erraticity behavior of the low multiplicity system, a 
new method based on measuring the rapidity gaps has been proposed by Hwa 
and Zhang [9]. They have proposed the moments 

n 1 n 

Gq = W+i ^ ’ ancl Hq = W+i ^ {1 ~ Xi) ~ q (3) 

' i=0 ' i—0 

as measures of spatial patterns in terms of rapidity gaps Xi = X i+ i — A, ; , i = 
0, • • •, N, with X 0 = 0 and X/v + i = 1. From this we can obtain the entropy¬ 
like quantities 


(G q In G q > 

< Gf In G f > ’ 


, v {H q \nH q > 
n q < H ( f In Hf > 


( 4 ) 


where G q and H q for experimental data, Gf* and for pure statistical 
fluctuations. How much degree of S q and Y> q deviate from 1 is a measure of 
erraticity in multiparticle production based on rapidity gaps. We calculate 
the S q and T, q only in central region (-2 < 77 < 2) and drop the events which 
has less than six particles in this region. A total of 2515 events for NA27 data 
and 2723 events for FRITIOF have been selected. 


The results for S q are shown in Fig.7a. It can be seen that S q deviate 
from 1 significantly and the InS'q versus q has a quite linear behavior. This 
means that S q satisfies the exponential relationship S q oc e aq . The straight 
lines are the linear fit to the experimental data. The fit parameter is listed in 
Table 3. This result is different from the power law behaviour claimed by Hwa 
and Zhang [9], S q cx q al , cf. the linear fit to the lnS^ versus lng plotted in 
Fig.7b and the fitting parameter listed in Table 3. The results from FRITIOF 
MC are also shown in Fig .8 (open circles) and Table 3. We can see that there 
is a same scaling behavior for the FRITIOF MC, although the values of S q 
deviate from experimental data sufficiently large. 


Table 3 The fit parameters obtained according to (13) and (14). 


EVENT SAMPLE 

a 

X 2 /NDF 

al 

x 2 /ndf 

experimental data 

0.133 ±0.002 

0.42 

0.66 ±0.06 

18.8 

FRITIOF 

0.078 ±0.002 

0.92 

0.38 ±0.04 

10.2 


The results for T, q are shown in Fig. 8 . It can be seen that T, q deviate from 
1 significantly, but lnS 9 versus q only has an approximately linear behavior 
S g oc e^ q , when q > 2. For the experimental data, (3 = 0.30±0.04. For the 
FRITIOF Monte-Carlo event sample, (3 = 0.15±0.03 . 

The value of f3 (0.15±0.03) for the FRITIOF MC is much smaller than 
that (/3 = 0.30±0.04) for the experimental data. The S q and T, q deviate 






















from 1 significantly implies that both of them are useful to serve as effective 
measures of erraticity in multiparticle production. 
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Figure Captions 

Fig.l. lniy versus ln(M rl M, f ) for different values of the Hurst-exponent. 
Fig.2 lnFjj versus InAf for continuous order q. 

Fig.3 D(q) versus q. 

Fig.4 f(a ) versus a. 

Fig.5 ln(l — /(«)) versus ln(R — a). 

Fig.6 The event space moments In C PA {M) versus In M. 

Fig.7 InS'q versus q and lng. 

Fig.8 lnE g versus q. 
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